Effects of free-electron-laser field fluctuations on the frequency response of driven 

atomic resonances 
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We study the effects of field fluctuations on the total yields of Auger electrons, obtained in the 
excitation of neutral atoms to a core-excited state by means of short-wavelength free-electron-laser 
pulses. Beginning with a self-contained analysis of the statistical properties of fluctuating free- 
electron-laser pulses, we analyse separately and in detail the cases of single and double Auger reso- 
nances, focusing on fundamental phenomena such as power broadening and ac Stark (Autler-Townes) 
splitting. In certain cases, field fluctuations are shown to influence dramatically the frequency re- 
sponse of the resonances, whereas in other cases the signal obtained may convey information about 
the bandwidth of the radiation as well as the dipole moment between Auger states. 

PACS numbers: 32.80.Aa, 32.80.Hd, 32.70. Jz, 32.80.Rm 



I. INTRODUCTION 

In traditional weak field photo-absorption processes, 
be it ionization, resonance excitation, scattering, etc., the 
only property of the source, besides the intensity, that en- 
ters the description is the bandwidth. The latter is due to 
random fluctuations present in any source, and depend- 
ing on the origin of the radiation, i.e. the nature of the 
process that gave rise to its production, it may be simply 
phase fluctuations or in addition intensity fluctuations. 
But in any case, for single-photon, weak field processes 
the total yield is sensitive to the bandwidth only, and 
possibly to the particular form of the source line-shape, 
depending on the type and details of the experiment. 

The situation changes dramatically, however, when the 
source is sufficiently intense to induce non-linear pro- 
cesses. Even the simplest case of 2-photon absorption 
introduces non-trivial departure from the above rule. 
Thus, the rate of 2-photon ionization, in lowest order 
perturbation theory (LOPT) is proportional, not simply 
to the square of the intensity, but rather to the second 
order intensity correlation function, while N-photon ion- 
ization is proportional to the Nth order correlation func- 
tion P, In yet another example, the strong driving of 
a single-photon resonant transition between two bound 
states, by anything other than an idealized monochro- 
matic sourcejbrings into play field correlation functions 
of all orders [3[ . 

The effect of field fluctuations on non-linear photo- 
interaction processes were discussed extensively during 
the late 1970s through the 80s, which led to the devel- 
opment of a number of theoretical approaches capable 
of addressing essentially any arrangement of atomic lev- 
els, including the continuum, coupled to one or more ra- 
diation fields involving fluctuations 3 8]. The relevant 
laser sources of that period, although pulsed, were of rel- 
atively long duration, in the range of nanoseconds to tens 
of picoseconds, as a result of which they could to a good 
approximation be viewed as stationary. Under the as- 



sumption of stationarity and crgodicity, certain standard 
models of stochastic fluctuations developed in quantum 
optics, such as a field with phase diffusion only, or a 
chaotic field (in the Glauber sense) Q , were employed in 
the development of the theories. Even so, the problems 
involved significant mathematical complexity, but even- 
tually, definitive, clear and in many cases rigorous results 
were obtained. In addition, a number of seemingly puz- 
zling experimental results were explained. 

With the advent of new sources of higher intensity, 
and much shorter pulse duration, in the few femtosec- 
ond range which in the infrared translates into even few 
cycles, the issue of field fluctuations ceased to be of rel- 
evance. That is because such short pulse sources are 
smooth, Fourier limited, whose bandwidth stems only 
from the duration, without any fluctuations of either 
phase or intensity [1 01 ] . Although frequency chirping may 
at times be present and even desirable, it can be reliably 
considered deterministic, without any stochastic aspects 
entering the description. 

The situation has changed again with the appearance 
of the short wavelength, accelerator-based free electron 
lasers (FELs), with photon energies ranging from the 
XUV to hard X-rays and pulse durations ranging from 
a few to a few hundred femtoseconds. Since the fre- 
quency in this regime is high, even a few femtoseconds 
correspond to many cycles of the field. Moreover, al- 
though the intensities in terms of W/cm 2 may be large, 
say 10 16 -10 18 , the photon flux and ponderomotive ener- 
gies are such that perturbation theory of the appropriate 
order is valid. This means that, depending on the situ- 
ation, rate equations or density matrix equations within 
a restricted atomic basis arc completely valid [ll| . Typ- 
ical cases in point are multiphoton ionization, which can 
also co-exist with a few-photon resonant excitation of a 
discrete state, known as REMPI (Resonantly Enhanced 
Multiphoton Ionization) [12, H3j |- But in this range of 
wavelengths, such discrete states will be unstable, de- 
caying by autoionization or Auger. Also, multiphoton 
ionization as well as resonant excitation at these large 
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photon energies, involve subvalence electrons; a feature 
that sets these processes apart from their counterpart in 
the optical regime, for which the above mentioned tech- 
niques, for treating field fluctuations, were developed. 

But with history repeating itself in some sense, the 
FEL pulses exhibit significant intensity fluctuations both 
in frequency and time domain; often referred to as spikes 
|14jla | . Given that typically the pulse, although short in 
terms of femtoseconds, it contains many cycles, and that 
it also involves amplified spontaneous emission, it could 
be argued that adopting the model of chaotic radiation 
should be an adequate approximation. In many cases and 
contexts, that may actually be satisfactory. But chaotic 
radiation, whose electric field classically is described by a 
complex Gaussian random variable, is a very well defined 
concept, which is strictly meaningful for a stationary 
source [9j. FEL radiation, however, cannot be assumed 
to be stationary, while its electric field, depending on the 
regime of operation, can be far from a Gaussian random 
variable. The theoretical handling of non-Gaussian ran- 
dom processes represents a formidable challenge. Leaving 
a more detailed description of the necessary models and 
tools for the next section of the paper, it will suffice to 
point out here that the analytical simplicity, available for 
a phase diffusion or a chaotic field, is no longer available, 
which implies a pure numerical simulation of the field 
and its stochastic properties. Such numerical simulation 
is not brute force, but based on detailed studies of the 
properties of the electron beam in the accelerator and the 
process of lasing taking into account the various regimes 
of operation of the machines. Our modeling reflects the 
state of knowledge of the FEL properties as of this time, 
in the self-amplified-spontaneous-emission (SASE) linear 
regime of operation. It is expected that with accumulat- 
ing experience and improving technology, control of the 
temporal properties of the pulses may improve, with the 
non-negligible possibility that in the future Fourier lim- 
ited pulses may be achievable. Until such time, however, 
the stochastic nature of the radiation needs to be taken 
into detailed account for its interaction with matter to 
be fully understood. On the other hand, it should be 
kept in mind that, the outcome of a non-linear interac- 
tion, depending as it does on the correlation functions 
of the field, provides a probe of the stochastic proper- 
ties of the source which are very difficult, if possible at 
all, to probe otherwise. In other words, the fluctuations 
inherent in the source, which for certain purposes may 
be a nuisance or dirt effect, from a different perspective 
contain information on the nature of the process that 
produced the radiation. That information, which goes 
well beyond the standard features of intensity and band- 
width, is revealed through the correlation functions. Put 
otherwise, a non-linear photo-interaction process probes 
the inner structure of the stochastic process underlying 
the generation of the radiation. 

Although as mentioned above, considerable theoreti- 
cal know-how, and rules of thumb, concerning the effects 
of field fluctuations on non-linear processes have accumu- 



lated over the years through studies in the optical regime, 
it is far from obvious that they can be applied intact 
when it comes to FEL radiation. In fact recent experi- 
ence in our own work connected to experimental results 
at FLASH, has led us to the reexamination of the funda- 
mentals of the theory, which has also been the motivation 
for this paper. As we show in the sections that follow, 
some of the accumulated wisdom is to a good approxi- 
mation transferable to this new regime, while significant 
modification and reformulation is necessary under cer- 
tain circumstances. And since our purpose here is to 
reexamine the fundamentals, we have chosen to focus on 
two basic themes. One is the simple single-photon ex- 
citation of an Auger resonance, the other being double 
Auger resonance in which one of the fields has the prop- 
erties of SASE FEL. The arrangement of DR, combined 
with the strong driving of either of the two transitions, 
with the other serving as a probe, appears in more than 
one context, as for example in EIT (Electromagnetically 
Induced Transparency) [l7[ or in the strong coupling of 
two highly excited resonances, embedded in a continuum, 
which provides an otherwise unavailable test of the the- 
oretical modelling of such resonances (l8l - [2l| . Our il- 
lustration of the effect of fluctuations on these arrange- 
ments through quantitative results and discussions, refer 
to examples taken from real experimental situations fl3| , 
which although specific, serve the purpose of calibration 
and extrapolation to essentially any analogous arrange- 
ment. 

Related studies in the past have not focused on the 
effect of fluctuations (22|, or have been concerned with 
the energy distribution (spectrum) of the Auger electrons 
[23l |24| . The aim of the present work is to study the ef- 
fect of fluctuations (inherently present in the pulses typi- 
cally produced in current SASE FEL sources) , on the fre- 
quency response of s ing le and double Auger resonances. 
By contrast to [23|, Hljj we do not analyse the actual 
spectrum of the Auger electrons, but rather we focus on 
the dependence of the total Auger yield on the detuning 
of the driving fields from resonance; a less demanding 
quantity that can be measured in related experiments. 
The following section is devoted to a brief description of 
our algorithm and its convergence. The single and dou- 
ble Auger resonances are discussed in Sees. UTT] and IIV1 
respectively. 



II. SIMULATION OF CHAOTIC SASE-FEL 
PULSES 

The statistical properties of light pulses emitted by 
a SASE-FEL depend crucially on the regime of opera- 
tion (l4j . In the regime of exponential growth (linear 
regime), it has been shown, that the pulses exhibit the 
properties of the so-called chaotic polarized light (13 - 
|16| . This means that the slowly- varying amplitude of 
the electric field is a complex Gaussian random variable, 
while the instantaneous intensity and the energy in the 
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pulses fluctuate according to known distributions. As the 
amplification process reaches saturation and the SASE- 
FEL enters the non-linear regime, the statistics of the 
radiation deviate significantly from Gaussian 14j . The 
theoretical description of non-Gaussian random processes 
is a rather difficult task, and to the best of our knowl- 
edge up to now no analytic results, for the properties of 
SASE-FEL radiation in the non-linear regime, have been 
obtained. Finally, a SASE-FEL operating deep in the 
non-linear regime seems to exhibit properties analogous 
to chaotic polarized light again [l4[ . 

Over the last decades, a number of FEL simulation 
codes have been developed, which take into account the 
details of particular FEL facilities [HI, HH and are thus 
capable of simulating its operation throughout the entire 
regime of parameters. Such codes, however, are rather in- 
volved and some times adapted to a specific FEL system. 
On the other hand, for the investigation of problems per- 
taining to the interaction of FEL radiation with matter, 
we need more direct ways to simulate the FEL radiation, 
or at least those properties that are pertinent to the par- 
ticular process under consideration. This is not always 
an easy task especially when the radiation involved does 
not obey Gaussian statistics. 

Throughout this manuscript, we consider light pulses 
produced from a SASE-FEL operating in the linear 
regime. The chaotic light is a fundamental concept of 
quantum optics, with the discussion usually limited to 
stationary and ergodic thermal sources (e.g., see p6j). 
In contrast to such a type of sources, SASE-FELs are 
not continuous sources, and they produce random light 
pulses which exhibit spikes both in time and in frequency 
domain fl4l - fl6l ] . Typically, the nominal duration of such 
a pulse is larger than the short time-scale of the field 
fluctuations (i.e., the coherence time). Such a type of 
radiation cannot be considered either ergodic or station- 
ary [27j . and thus simplifications and analytic expres- 
sions typically used for thermal sources do not apply in 
the present scenario. For instance, ensemble averages 
of time-dependent quantities cannot be substituted by 
integrations over time. In the following we adopt numer- 
ical techniques that have been developed in the context 
of quantum optics p8l - l30| , in order to produce fluctuat- 
ing pulses, which exhibit all the properties of SASE-FEL 
pulses in the linear reg ime. The details of our algorithm 
can be found in (2!|[30[ but, for the sake of completeness, 
we briefly summarize here its main aspects. It should 
also be noted that our algorithm bears analogies to the 
algorithm used by other authors (23|. 



A. Algorithm 

In various contexts of physics, one usually deals with 
noises of different origin, i.e., with various types of fluc- 
tuating random variables. Gaussian white noise typically 
refers to a quantity that fluctuates e.g., with time, and it 
has flat power spectral density (PSD). Formally speaking, 



it is represented by a random variable with Gaussian 
distributed fluctuations, which satisfies 



(£(*)) = <e 2 w) 



= o, 
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where S(-) is the Dirac delta function, while (•) denotes 
an ensemble (statistical) average. Such a noise can be 
generated e.g., from uniform random variatcs by means 
of Box-Miillcr algorithm. 

The generation of Gaussian colored noise, however, is 
far from straightforward and has been the subject of a 
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rather limited number of papers (e.g., see 
all of these papers, the algorithms take as an input uni- 
form or Gaussian white noise, and return Gaussian noise 
£(i) with autocorrelation different from the Dirac delta 
function i.e., the new variable £(t) satisfies Eq. (|lal) . 

but (C(i)C(i')) ~ with G(\tit') ^ S(t - t'). 

The colored- noise algorithm of [28|, relies on the propa- 
gation of a stochastic differential equation, and is limited 
to the generation of exponentially correlated Gaussian 
colored noise (Ornstein-Uhlenbeck process). Through- 
out this work, we adopt the algorithm of [29|, |3(| , which 
is capable of producing Gaussian colored noise for various 
types of autocorrelation functions, from the correspond- 
ing PSD, utilizing the fast Fourier transform. The perfor- 
mance of the algorithm has been discussed in the context 
of particular stochastic differential equations [29|,[30|], and 
here we will discuss briefly the performance in the con- 
text of simulations pertaining to SASE-FEL pulses. 

The algorithm is implemented on a grid of N g points 
in frequency domain with step 5u), from uj s — A/ 2 to 
uj s + A/2, with A = NgSuj and oj s denoting the central 
frequency of the PSD &q(uj) under consideration. An 
independent complex Gaussian random variable cor- 
responds to the fcth point of the grid (of frequency u>k) 
, with (fr) = (e k ) = and = 5uPt{u k )6 k , u 

where 8^,1 is the Kronecker's delta. The generation 
of noise in the time domain is achieved by means of 
the discrete Fourier transform, obtaining complex Gaus- 
sian random variables with autocorrelation function 

Gf\t,t') = {((%*(?)), given by 

N g /2 

Gf\t,t')= Sue^-^&cfa). (2a) 

k=~Ng/2 

It has to be emphasized here that this algorithm gener- 
ates a stationary random process described by the ran- 
dom variable ((t), with (£(t)) = and 
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(2b) 



In the context of SASE-FEL, one is interested in the gen- 
eration of fluctuating pulses i.e., signals of finite duration 
and of time-dependent average profile. To this end, the 
noise produced by means of the aforementioned proce- 
dure can be superimposed to a particular profile (enve- 
lope) yielding thus a spiky pulse, analogous to the ones 
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FIG. 1. (Color online) (a,c) The dotted curves show two of 
the profiles f s (t) used in our simulations. (b,d) A sample 
of two random spiky pulses, typically produced in a single 
realization of the algorithm discussed in Sec. [Ill by super- 
imposing Gaussian correlated noise OJ 1 = 7.2 fs) with the 
deterministic profiles of (a,c). The solid curves in (a,c) show 
the average intensity (I s (t)) //J on a sample of 1000 random 
spiky pulses. 



typically produced from SASE FEL sources. Clearly, the 
new signal cannot be considered stationary. 



B. Application and Convergence 

The above algorithm has been tested for various types 
of colored noises and various profiles. In the following 
we present some results pertaining to the generation of 
fluctuating pulses with Gaussian correlated noise, which 
is the type of correlations typically observed in different 
SASE FEL facilities (TlM M, M, HI . Since our noise 
has to be Gaussian correlated, the algorithm of Sec. Ill A[ 
is seeded with a Gaussian PSD for the noise i.e., 



exp 



(u - u s y 



2a, 2 



(3) 



where a u is the standard deviation of the distribution. 
As will be discussed later on, the generated noise £(i), is 
a complex Gaussian random variable (i.e., the real and 
imaginary parts have Gaussian distributions), and ex- 
hibits Gaussian correlations, with the coherence time T c 
determined by er w . 

To facilitate our theoretical treatment, the amplitude 
of the electric field with central frequency to a in a single 
realization of the algorithm (within some non-essential 
multiplicative constants) is defined as 



(4) 



where f s (t) is a Fourier- limited (deterministic) pulse 
profile of finite duration and peak value J s . The inten- 
sity of the stochastic pulse in the time domain is simply 



I s (t) = \£ s (t)\ 2 =I^f s (t)\C(t)f 



(5) 



The deterministic envelope f s (t) ensures the smooth rise 
and drop of the intensity, and can be chosen at will. In 
Fig. Q]wc show two of the profiles used for f s (t) through- 
out our simulations, together with a small sample of spiky 
pulses. The profile of Fig. QJa) is a Gaussian given by 



f s (t) = exp 



(6) 



where r s is the pulse duration and to > is the center of 
the pulse. The profile of Fig. [IJc) is a mathematical con- 
struction obtained by superimposing different Gaussians, 
and enabled us to check the dependence of our results on 
the shape of the profile. From now on we refer to this 
profile, as profile 2. In view of Eqs. (|2bp , (|3|) and (|5|), 
averaging over a large number of random pulses (which 
means random functions), one recovers the deterministic 
profile f s (t) i.e., 



(I s (t))=li 0) f s (t). 



(7) 



Typically, one has to average over at least 10 3 random 
spiky pulses to ensure a rather good convergence in this 
respect [sec Fig. [TJa,c)]. Throughout our simulations 
time evolution takes place for t € [0,1/], with 7) a fi- 
nite positive number and the pulse parameters (center 
and duration) chosen such that the pulse rises and falls 
smoothly in this time interval (i.e., formally speaking 
I s (t) -> as t -> {0,7)}). 



1. Statistics and correlations 

The fluctuations of the instantaneous electric field (not 
shown here) obey a Gaussian distribution. The proba- 
bility distribution of the instantaneous intensity, i.e., the 
probability for an instantaneous measurement to yield a 
value between I s (t) and I s (t) + dl s (t), is described by the 
negative exponential probability density function (PDF) 



rr / M 1 f hit) 



(8) 



The energy in a random pulse at a space point in the 
interaction volume is given by 



W s cx 



/«(t)dt, 



(9) 



and it fluctuates from pulse to pulse. The corresponding 
probability distribution is the Gamma PDF given by 



p{W s ) 



M M W M-l 

'—M eX P 



r(M) (Ws 



-M 



Ws 



(10) 



where T(M) here is the Gamma function with argument 
M = (W s ) 2 I {(W s - (W s )) 2 ). The parameter M repre- 
sents the average number of modes in a radiation pulse. 
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All of the above properties are in agreement with experi- 
mental observations and theoretical results pertaining to 
various SASE-FEL facilities [TIMS!- Although, ow- 
ing to the inevitable focusing in experiments with high 
intensity sources, the interaction volume relevant to a 
particular experiment needs to be taken into account in 
the quantitative interpretation of data, this is a separate 
issue. The temporal characteristics of the radiation are, 
however, the same throughout the interaction volume, 
which means that our results can be applied to any ex- 
perimental arrangement by folding in the geometry of the 
interaction volume (see also related remarks in Sec. fV|) . 

The negative exponential distribution for the inten- 
sity and the Gamma distribution for the energy are inti- 
mately connected to a field that obeys Gaussian statistics 
[2(|[27j]. The corresponding first-order coherence, how- 
ever, is not related to the Gaussian statistics of the field, 
but rather to the chosen PSD of the generated noise. 
Choosing a Gaussian PSD of the form (J3|), one expects 
a Gaussian first-order autocorrelation function. By defi- 
nition G {1) (t,t') = (£ s (*)£*(*')), and using Eqs. (0} and 
([7]) we obtain 

|G«(M')I = y/UM) (W))\G^(t,t')\ 

= 4°Wf s (t)fs(t')\Gf ) (t,t')\, (11) 



where (*,*') 



f\t,t') = (C(t)C* (*'))• We see therefore that the 
first-order autocorrelation function for the field not only 
depends on the statistical properties of the noise ((t), but 
also on the ensemble average of SASE-FEL pulses (I s (t)) . 
Only in the case of stationary fields one has constant 
(I s (t)), and thus G^\t,t') depends solely on the noise 
correlations and not on the average profiles (we return to 
this point at the end of the section). The modulus of the 
degree of first-order temporal coherence is defined as [26| 



|9 (1) (M')I 



^(I s (t)) (l a (t')Y 



(12) 



and is equal to \Gf\t,t')\. For the Gaussian PSD (©, 



we have 



\g^(t,t')\ =cxp 



alit-t'Y 



= \G?>(t,t% (13) 



with cr w determining how fast the correlations drop with 
the delay. 

For a field that obeys Gaussian statistics, Wick's the- 
orem implies that higher-order degrees of coherence can 
be expressed in terms of g^(t,t'), and hence they are 
determined by the spectral properties of the field (26j . 
Throughout our simulations we have used various types 
of ^c(w) and thus of correlations, some of which are 
summarized in table U together with the corresponding 
bandwidths 7 [i.e., the FWHM of ^(w)]. The coherence 
time is defined as 



7«f 



\9-'(v)\ 2 dv, (14) 

where v = t — t' , and is also shown in table HI For the 
discussion to follow, it is worth keeping in mind that for 
fixed <j u and for v < 2, the exponential I^^H^)! exhibits 
the fastest drop, and the hyperbolic secant the slowest 
one. 



2. Bandwidth of fluctuating pulses 

The bandwidth 7 is determined by the noise in the 
pulse, and a question arises here as to the bandwidth of 
the fluctuating pulses produced by superimposing this 
random noise with the Fourier-limited average profile 
(I s {t)). This profile ensures the smooth rise and drop 
of the intensity for t G [0, 00) and thus we can safely as- 
sume that I s (t) is a square integrable function. One can 
then define the energy spectral density of the random 
pulses as S s (u>) = (|£ s (w)| 2 ), where £ s (u>) is the Fourier 
transform of £ s (t), and using Eqs. Q and (JT)) we obtain 



(2tt)2 J_ 



00 poo 



dtdt'y/(I a (t)) (I a (t'))G^\t,t')e 



(15) 



00 j — 00 



This is the Wiener-Khintchine theorem for non- 
stationary square integrable signals. For given (I s (t)) 
and specified noise, one can obtain (a s (oj), and thus the 
associated bandwidth. Usually one is interested in the 
normalized spectrum S s {uS), which is obtained by divid- 
ing S s (w) by J^S's^duj. 

For Gaussian (I s (t)} of duration r s , and the Gaussian 
correlated noise of Eq. (|13|) , <S S (ui) is also Gaussian. We 
obtain 



- , . 2Vln(2) r 41n(2)w 2 l . . 

£ a {w) = f-^ ; exp ±±- , (16a 



where Auj s is the bandwidth (FWHM) of & s (u>) given by 
Aids = 2 V /TT ^V /I ^), (16b ) 

Ts 

The case of Fourier-limited pulses corresponds to 
<t u t s -c 1, with the result 

A.- = = (16c) 

t s At s v ' 

which is the time-bandwidth relation for Gaussian 
Fourier-limited pulses usually found in standard text- 
books with the FWHM of (I s (t)) denoted by At s = 
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TABLE I. Various PSDs with the corresponding degrees of first order temporal coherence, bandwidths, and coherence times. 



Label 


£» c (w) a 


[\9 W (v)\t 


7 


T c 


Exp. 




exp[— v] 






Gauss. 


cxp[-w 2 /2]/(^\/2^) 


exp[-C 2 / 2 ] 


2a uV /21n(2) 




Sech. 


sech(7riD) /cr u 


sech(5/2) 


20^ cosh -1 (2) /7T 


^j 1 



b v = (t - t')cr, 



2y/\n{2)r s . Hence, Eq. (|16bl) may be also expressed as 
&lo s = AwfVl + V, (16d) 

where 

X = OujT s - (16e) 

In the other limit, even for x > 2, one essentially ap- 
proaches the limit of stationary source with Aoj s ~ 7, 
where 7 is given in table HI As mentioned above, a u is 
inversely proportional to the coherence time of the noise 
entering the pulse (e.g., see table IJ), and thus % charac- 
terizes the number of spikes in a Gaussian pulse profile 
of duration r s . The derivation of simple analytic expres- 
sions for <t> s {uS) is a rather difficult task, when the noise is 
not Gaussian correlated and / or for arbitrary average pro- 
files. In any case, however, the bandwidth of a fluctuating 
pulse Aw s is expected to be determined by the bandwidth 
of the Fourier-limited average pulse profile (I s (t)) and the 
product x, albeit perhaps through a rather complicated 
expression. From another point of view, expressing a u in 
terms of 7 (see table H| , Eq. (|16dj) reads 

A. s = A.ry / l + (^j" • (17) 

This relation shows that for pulses with Gaussian average 
profile that exhibit Gaussian-correlated fluctuations, the 
bandwidth is the geometric mean of the bandwidths cor- 
responding to the Fourier-limited average profile and the 
fluctuations. For 7 ^> Aw™°, the bandwidth of the pulse 
is fully determined by the fluctuations i.e., Aw s ~ 7, 
whereas for 7 -c Aw™, we have the case of a Fourier- 
limited pulse with Aw s ~ Aw™ n . 

3. Convergence with respect to intensity moments 

Let us discuss briefly the convergence of the algorithm 
in terms of the moments of the instantaneous intensity. 
In view of Eq. (O we have, 

/>oo 

(J.(t) r ) = / I s (t) r p[I s (t)Ws(t) = r! (I s (t)) r . (18) 
Jo 

As depicted in Fig. [2j the ratio (I s (t) r ) / (I s {t)) r fluctu- 
ates around H, with decreasing amplitude of oscillations 
as we average over more pulses (realizations) [Hj]. For 
increasing r, the convergence is slower, in the sense that 




FIG. 2. (Color online) Convergence of the algorithm discussed 
in Sec. [Ill The ratio (Is(t) r ) / (I s (t)} r is plotted as a function 
of time, with the average taken over two different numbers of 
pulses (realizations). The thick grey (horizontal) lines refer 
to r!, which is the theoretically expected value of the ratio. 
Other parameters: r s = 10 fs, crj 1 = 2 fs. 



more realizations are needed to approximate adequately 
rl. Note, however, that even for 10 3 or 5 x 10 3 real- 
izations, the amplitude of the oscillations for r < 5 is 
smaller than the gap between successive factorials, and 
thus one can distinguish unambiguously between differ- 
ent values of r in this regime. The appearance of the r! 
factor in Eq. (fT5| is a distinctive feature of chaotic light 
[26|, which is expected to affect considerably fundamen- 
tal nonlinear physical processes, such as the multiphoton 
atomic ionization II il l. 



4- Stationarity 

The aforementioned algorithm can be used in stud- 
ies pertaining to interactions between matter and SASE- 
FEL radiation with the statistical properties described 
above, and various types of autocorrelation functions. In 
the following we explore the influence of fluctuations on 
the frequency response of atoms. Analogous studies in 
the past have focused on stationary fields i.e., for time- 
independent (I s (t)) [1, 0,13 1 while most of them consid- 
ered exponentially-correlated fluctuations. In the present 
SASE-FEL facilities such conditions are not satisfied, and 
the ensemble average intensity (I s (t)) has a finite dura- 
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FIG. 3. (Color online) An atomic transition pertaining to 
a sub-valence electron, driven by chaotic SASE-FEL pulses. 
The resulting inner-shell vacancy in the excited atom decays 
with rate IV 



tion t s . It may be possible to assume some sort of station- 
arity only if the coherence time T c is much smaller than 
the pulse duration t s , and the field is observed for a time 
window much shorter than r s . Such a condition may be 
fulfilled in practise when one monitors directly the light, 
but it is hard to be fulfilled when matter interacts with 
SASE-FEL pulses. Typically, in such cases the target 
(atoms, molecules, etc) experience the rise as well as the 
fall of each one of the pulses, and one simply monitors 
the products of the interaction (i.e., electrons, ions, etc). 
Under such circumstances, the SASE-FEL radiation may 
or may not be considered stationary and any theoretical 
description has to take into account the finite temporal 
width of (I s (t)). As long as {I s {t)} is a smooth func- 
tion, its details besides the duration r s , are not expected 
to play a significant role in most cases. For the prob- 
lems under consideration, we have confirmed this fact by 
considering different profiles for (I s (t)), two of which are 
shown in Figs. HJa) and (c). 



III. SINGLE AUGER RESONANCE 

Consider the case of a sub- valence atomic transition 
depicted in Fig. [3] A SASE-FEL beam is focused on a 
target of neutral atoms, inducing an clcctric-dipole tran- 
sition of an inner-shell electron from state |1) to a highly 
excited state |2). The relaxation of the sub- valence va- 
cancy in the excited atom via an Auger decay, gives rise 
to Auger electrons which are observed in the experiment. 
Let hujj denote the energy of the state and let T 2 be 
the rate associated with the Auger decay. The electric 
field of the radiation is given by 



E a (t) = [s s (ty 



(19) 



where lu s denotes the central frequency of the spectrum, 
e s is the polarization vector, and £ s (t) is the fluctuating 
complex amplitude. 



Throughout this work £ s (t) is treated as a stochas- 
tic complex Gaussian random function and is generated 
along the lines of the previous section. The instantaneous 
Rabi frequency O s (f), given by 



fl 12 ■ e s £ s (t) 



(20) 



is also a stochastic complex Gaussian random function 
with zero mean and variance determined by the variance 
of the field. In this definition, fl is the electric dipole 
operator, and ff\2 is the transition dipole moment for 
|1) <-> 1 2). In the following for the sake of brevity we also 
write /ii2 = M12 • e s . In view of Eqs. ([4]) and ([7]), we have 



Pl2 



C(t) = n^^WK(t). (2i) 



where 



Pl2< 



(22) 



is the peak value of the Rabi frequency. 

The problem can be formulated in the framework of 
the reduced atomic density matrix with elements pij (t) . 
In the rotating- wave approximation, the equations of mo- 
tion for read 



dan 

dt 
da 2 2 
~dT 
d<J\2 

dt 



2Im[0> 12 ] 

-r 2 (722 - 21m [0*012] 



(23a) 
(23b) 



iA s - ) (T12 + ifi s ((T22 - on), (23c) 



where p u = an, p 12 = cr^e 1 ' 



A., 



^21 



is the 



detuning of the field from resonance, while Stark shifts 
have been neglected. In the absence of other types of 
(in) homogeneous broadening mechanisms, we have Ti 2 = 
T2. In the presence of fluctuations in the electric field, 
this is a set of coupled stochastic differential equations. 

The frequency response of the atoms to the SASE-FEL 
radiation, as we vary u> s around resonance, is obtained by 
monitoring the Auger electrons. The total probability for 
Auger decay is given by 



= r 5 



0-22 (*) dt. 



(23d) 







Alternatively, we can add to the equations of motion for 
<7ij , the following differential equation 



dQ 2 
dt 



r 2 o 2 2, 



(23e) 



where T 2 is the probability per unit time for Auger decay. 

In each realization of our simulations a SASE FEL 
pulse £ s (t) is generated randomly according to the al- 
gorithm described in Sec. [TTJ Subsequently, Eqs. (|2"3"|) arc 
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FIG. 4. (Color online) Single resonance driven by a Fourier- 
limited pulse. The Auger signal Q2 is plotted as a function of 
the detuning A s . (a) Various values of the ratio fii ' /T 2 and 
for t s — SF^ 1 ; (b) Various values of T2T S for fixed = 2F 2 . 
A close-up of (a) and (b) around the peak is shown in (c) and 
(d), respectively. The dashed curves show the best Lorentzian 
fits to the numerical data for (b) F2T S = 4 and F2T S = 8; (c) 
fli ^ — 0.5r2 and fl s ^ — 4.0r2. Other parameters: Gaussian 
pulse profile, T) = 32F2" 1 . 
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FIG. 5. (Color online) Single resonance driven by stochastic 
pulses. Stochastic Auger signal averaged over 5000 stochastic 
pulses as a function of the detuning A S) for various values of 
the ratio fl s 0) /r 2 and (a) X = 1-67; (b) X = 2.5; (c) X = 5.0; 
(d) X = 10.0. For the sake of comparison, the signal for a 
Fourier-limited Gaussian pulse with intensity profile the same 
as (I s (t)) is also shown for = 0.5r2 and f)i°' = 2F2 

(dashed curves). Other parameters: Gaussian profile {I s (t)) 
with t b = Sr^" 1 , Tf = 32F^" 1 , Gaussian-correlated noise. 



propagated from t = to t = Tf, keeping track of Q2 at 
t = Tf, where Tf ~S> t s ; which is essentially equivalent to 
taking the upper limit of Eq. (|23d[) to infinity. Averag- 
ing over a large number of random pulses, we obtain the 
average stochastic signal (Q2) ■ 

In connection to the scheme of Fig. [3l it is worth men- 
tioning here the recent experiment by Mazza et al. [Hj], 
pertaining to the spectral response of the Auger reso- 
nance 3d— >5p in Kr. This experiment was in the regime 
of weak fields and long pulses (with respect to T^ 1 ). By 
contrast, the pulse durations considered throughout the 
present work are comparable to T^ 1 , while a broad range 
of peak values for the Rabi frequency is considered. To 
keep our formalism as general as possible, in the following 
discussion the various quantities are measured in units 
of the natural linewidth T2 . The actual value of T2 de- 
pends on the particular atomic resonance under consid- 
erations. For instance, the linewidth for the Auger res- 
onance 3d— >5p in Kr is 83 meV, which corresponds to a 
lifetime of about 8 fs [ljj. The Rabi frequencies we con- 
sider here for the FEL radiation vary from 0.1 to 2T2, 
corresponding to intensities 10 14 — 10 17 W/cm 2 for the 
3d— >5p resonance in Kr, whereas the intensities in the 
experiment of [l3| were much lower. 



B. Numerical results 

The behaviour of Q2(A S ) for a Gaussian Fourier- 
limited pulse of various durations and various peak Rabi 
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FIG. 6. (Color online) Single resonance driven by stochastic 
pulses. Stochastic Auger signal averaged over 1000 stochastic 
pulses as a function of the detuning A s for (a) Gaussian- 
correlated noise and various combinations of pulse durations 
and coherence times; (b) various types of fluctuations and 
t s = SF^ 1 , = 10T7 1 . The curves of (a) correspond to the 



combinations (from 1 to 4): 
-1 



Ts — 3r 2 and cr w = 10r s ; t s = 



4.5r 2 and a u 



3F„ and Uu. 



20r; 



4.5F2 -1 and a u = 20r s _i , respectively. The degrees of first- 
order coherence for the fluctuations chosen in (b) are shown 
in table [TJ For the sake of comparison, in Fig. [6] we also plot 
the signal corresponding to a phase-diffusion model (PDM) 
0) HI with 7 = 2^/21n(2)a„. Other parameters: Gaussian 
profile for (I s {i)), n s 0) = 2r 2 , T f = 32F2 1 . 



frequencies is depicted in Fig. |U As depicted in Fig. 
|Ua), for fixed t s and T2, Q2(A S ) becomes wider with 
increasing peak Rabi frequency f2s . The derivation of 
an analytic expression for Q2 is a rather difficult task be- 
cause of the time dependence of fl s (t). Nevertheless, our 

simulations show that for weak fields (small O £° ^ ) , the de- 
pendence Q2(A S ) is well approximated by a Lorentzian. 
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This is the case of the experimental setup in [13|. For 
increasing f2i°\ however, we observe here severe devi- 
ations from the Lorcntzian, especially around the peak 
of Q2(A S ) [see Figs. 0Ja,c)]. Furthermore, according to 
Figs. |U[b,d), Q S (A S ) becomes more flat at the apex, 

as we increase the pulse duration t s , for fixed fti°^ and 
r 2 . The atoms experience large intensity for longer pe- 
riods as we increase r s , and thus the deviations from 
the Lorentzian extend to even larger dctunings, whereas 
the Lorentzian seems to remain a good approximation 
for large A s only. These observations do not depend on 
the actual pulse profile under consideration, and can be 
attributed to power broadening under the constraint of 
Q2(A S ) < 1 for all values of A s . Note, however, that 
in view of the pulsed driving, the power broadening de- 
pends on the peak intensity as well as the duration of the 
pulse. 

In Fig. [5J we show results for the average stochastic 
signal (Q2) as a function of A s , in the presence of fluctu- 
ations. Recall here that the bandwidth of the stochastic 
pulses with Gaussian average intensity profile exceeds the 
bandwidth of a Fourier limited pulse of the same profile 
and duration by ^/l + 2\ 2 - The duration of the aver- 
age intensity profile has been chosen t s = Sr^" 1 and the 

peak Rabi frequency flf^ < 2r 2 . According to Fig. HJ 
the deviations from the Lorentzian shape are negligible 
for a Fourier-limited pulse of the same parameters. As 
depicted in Fig. [3J however, the presence of amplitude 
fluctuations in the pulse enhance the power broadening 
and thus, strong deviations from the Lorentzian appear 
even for the particular combinations of t s and Qi . For 
X ~ 1 the coherence time of the noise is comparable to 
the duration of the pulse [see Fig. |5ja)], and (Q S (A S )) is 
very close to the signal for a Fourier-limited pulse with 
the same profile as (I(t)). Increasing x, we essentially de- 
crease the coherence time relative to r s thus approaching 
the limit of stationary source (Alu s w 7). Accordingly, 
the average stochastic signal (Q 2 ) becomes significantly 
broader than the signal for a Fourier-limited pulse of the 
same parameters, and only for large values of A s , the 
two signals approach each other [compare the dashed thin 
curves to the corresponding thick coloured curves in Figs. 
[SJb-d)]. Note also the prominent flattening of (Q 2 ) at the 
apex for x ^> 1 [sec Figs. [5jc,d)]. 

Figure [S{a) shows the average stochastic signal as a 
function of A s , for various combinations of r s and a^. 
Increasing the duration t s , while keeping x constant, re- 
sults in a faster drop of (<3 2 (A S )) as we move away from 
A s = 0, and thus to a narrower Q2(A S ) (compare curve 
1 to curve 2, and curve 3 to curve 4). On the other hand, 
if we increase x f° r fixed r s , (compare curve 1 to curve 
3, and curve 2 to curve 4), the effect is the opposite. 
In Fig. E]L>), we plot (Q 2 (A S )) for fields with different 
types of noises (the corresponding PSDs and g^ (t — t') 
can be found in table |T| , but with the same . Clearly, 
in all cases (Q2) flattens in the neighbourhood of A s = 0, 
and the type of correlations affects mainly the wings of 
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FIG. 7. (Color online) Double Auger resonance driven by 
a short- wavelength SASE-FEL pulse (lower transition), and 
a Fourier-limited pulse (upper transition). The sub- valence 
vacancy in the excited atom decays via two distinct Auger 
channels with rates F2 and r$. 



(Q2{A S )), as we move away from A s = 0. The faster 
g^'fy—t') drops, the slower (Q2(A S )) drops with increas- 
ing A s . Moreover, for a field with phase fluctuations only, 
(Q 2 (A,,)) drops much slower than for a stochastic field 
of the same bandwidth with both amplitude and phase 
fluctuations (compare PDM to Gauss, curves). 

The above discussion and findings hold for all the dif- 
ferent profiles we used for (I(t)) throughout our simu- 
lations, albeit perhaps with some qualitative differences. 
We emphasize once more that the above discussion on the 
power broadening and the non-Lorentzian signals pertain 
to the total Auger yield as a function of the detuning of 
the driving FEL radiation from resonance. To the best 
of our knowledge, no related experimental observations 
have been reported in the literature so far. Perhaps the 
reason is that most of the related experiments have been 
performed in the regime of weak intensities. Our theory 
suggests that the required Rabi frequencies should be at 
least twice the natural linewidth T 2 of the observed Auger 
resonance. For the 3d— >5p resonance of Kr, this corre- 
sponds to intensities of the order of 10 17 W/cm 2 . Such 
intensities are certainly within reach of various FEL fa- 
cilities such as FLASH and LCLS. 



IV. DOUBLE AUGER RESONANCE 

Consider now the case depicted in Fig. [JJ where the 
core-excited state |2) not only decays, but it is also cou- 
pled to another state |3), via second field Ed(t). 
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In contrast to E s (t), we assume that Ed(£) is a Fourier- 
limited pulse (typically of a longer wavelength) given by 



E d (i) = [f d (t)e 



■£* d (t)e 



-ioJdt~] 



sd 3 



(24) 



where aj d denotes the frequency and e s is the polariza- 
tion vector. The complex amplitude £d(t) can be written 

as £d = \J I d °^ fd(t), where fd(t) is a Fourier- limited en- 
velope of peak value 1, I d °^ is the the peak intensity of 
the pulse, and let T d denote the corresponding duration. 
The transition |2) O |3) is assumed to be electric-dipole 
allowed, and in analogy to the definitions of the previous 
section, the corresponding Rabi frequency is given by 



n d (t) 



P2 3 sjld(t) _ (0) 



where 



n 



(0) 



^231 



(25) 



(26) 



is the peak Rabi frequency, and jl 23 the relevant transi- 
tion dipolc moment. The corresponding detuning from 
resonance is A^ = uj 32 — w<j, while we must also include 
the Auger decay from state |3) with rate r 3 . This is 
the so-called double-resonance (DR) setup which allows 
for the observation of ac Stark (Autler-Townes) splitting, 
and it is also employed in various quantum optics phe- 
nomena, such as the EIT [l7j |. 

The interaction of the atomic system with the radi- 
ation can be described by means of the reduced atomic 
density matrix with elements Pij(t). In the rotating- wave 
approximation, the equations of motion for pij read 



= 2Im[fi*<x 12 ] 

(9(799 

-T 2 a 22 - 21m [fi*<7i 2 - fl* d a 23 ] 



dt 

d<?33 
dt 

dai 2 
dt 

dcr 23 
dt 

d(J\3 
dt 



(27a) 
(27b) 

= -r 3 a 33 - 21m [ft> 23 ] (27c) 
= (iA s - T 12 )a 12 + 1 [Q s (a 22 - <j n ) - Q>i$7d) 
= (iA d - T 23 )a 23 + i [n*a 13 + n d (<r 33 - V22§7e) 

= [i(A s + A d ) - r 13 )]<T 13 + i (n s <j 23 - n d CTi#7f) 



where p u = a u , p 12 = a 12 e lu -\ p 23 = a 23 e iwit , 
P13 = ai 3 e l ^ 3+Ud ^ t . In the absence of other types 
of (in)homogeneous broadening mechanisms the off- 
diagonal relaxation rates are = (1^ +Tj)/2. 

The frequency response of the atoms is obtained by 
monitoring the Auger electrons. The total probability 
for Auger decay from state \j) is given by 



'.1 j 



(t) dt. 



(27g) 



In each realization of our simulations a SASE-FEL pulse 
£ s (t) is generated randomly, according to the algorithm 
described in Sec. HJ and the above set of stochastic differ- 
ential equations is propagated from t = to t = Tf with 
Tf 3> max{ T Sl r d }, keeping track of Qj at t = Tf. The 
average stochastic signal (Qj) is obtained by averaging 
over a large number of pulses. 

One can distinguish between two different arrange- 
ments of weak and strong pulses in the context of DR. In 
the first arrangement, the Fourier- limited pulse is strong 
(pump) and the stochastic SASE-FEL pulse is weak, 
and is used as a probe for monitoring the upper tran- 
sition. This arrangement is pertinent to EIT and has 
been discussed recently for optical control of X-ray ab- 
sorption [l7| . Alternatively, one may also consider a 
strong stochastic pump in the lower transition, and a 
weak Fourier-limited probe in the upper one. We have ex- 
plored both arrangements and in the following we present 
our main findings on the effect of fluctuations in the DR 
setup. To facilitate the characterization of such effects, 
in the appendix we summarize briefly some of the funda- 
mental aspects of the problem in the absence of fluctua- 
tions (Fourier-limited pulses). 



B. Numerical results 

In this section we discuss the effects of fluctuations in 
the DR signal. The combinations chosen for peak Rabi 
frequencies and durations for the pump/probe pulse pro- 
files are such that, in the case of Fourier-limited pulses, 
Qj exhibits a clear doublet when plotted as a function 
of the detuning of the probe from resonance. The deter- 
ministic signal Qj for Fourier-limited pulses serves as a 
reference for the investigation of effects induced by the 
fluctuations in the probe and the pump. In other words, 
the average stochastic signal (Qj) in the presence of fluc- 
tuations, is compared to the deterministic signal Qj per- 
taining to Fourier-limited pump/probe pulses. 



1. Arrangement I: stochastic probe 

Although usually one would want the probe not to have 
fluctuations, in considering FEL as a probe we are com- 
pelled to consider the effect of fluctuations (e.g., this is 
the case of the experiment [l7|). The stochastic probe 
is applied on |1) <R- |2), and a Fourier- limited pump on 
|2) o |3). The average stochastic signal (Q 2 ), is plotted 
as a function of A s in Fig. [51 for various values of x an d 
two different dctunings of the pump A^. For the sake of 
comparison, the deterministic signal for Fourier-limited 
probe pulse with profile (I s (t)) is also shown (dashed 
curve). For \ ~ 1, the signal is very close to the one 
for Fourier- limited probe. As we increase x (i- e -i as we 
decrease T c ~ <j~ x relative to r s ), the two peaks of the 
doublet become broader, whereas their position is also af- 
fected slightly. Recall here that according to Eq. (|16dl) . 
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FIG. 8. (Color online) DR arrangement I with stochastic 
probe and Fourier-limited pump. The average stochastic sig- 
nal is plotted as a function of the probe detuning A s for in- 
creasing values of x> an d for (a) resonant (Ad = 0) and (b) 
off-resonant (Ad = 2T2) pump. The thin dashed curves corre- 
spond to a Fourier limited probe (i.e., x <C !)■ Synchronized 
Gaussian profiles have been chosen for the average probe in- 
tensity (Is(t)) as well as for the pump. Other parameters: 

r a = 4.5r 2 ~ \ r d = 6r 2 -\ n< 0) = o.ir 2 , = ior 2 , r 3 = r 2 , 

Tt = 32rr 1 , 5000 realizations, Gaussian-correlated noise. 



Aw s w 7 for x > 2, with the bandwidth 7 given in the 
second row of table |U Thus, the visibility of the dou- 
blet decreases with increasing bandwidth (or decreasing 
T c ), and for x > 20 it has essentially disappeared. It 
is worth noting here that this critical value of x above 
which no doublet is observed corresponds to 7 « f^. In 
the picture of eigenstates |±) (see appendix), the fluctu- 
ating probe leads to fluctuating couplings |(1|±)|, which 
in turn give rise to a noisy signal, i.e., to broader peaks. 
When this broadening exceeds the peak Rabi frequency 
of the pump, the doublet disappears. 

In Figs. Eta-c), we plot the separation of the peaks 
in the doublet as a function of x" 1 ! f° r various types 
of noises (see table |T| . For the sake of comparison, the 
separation is normalized to its value for Fourier-limited 
probe with the same profile as (I(t)) (see dashed curve 
in Fig. [5J. As was expected, in all cases the normal- 
ized splitting approaches 1 as x — > 1, and decreases as 
we decrease x _1 (i-e., as we decrease T c ~ a~ x for fixed 
t s ). The decrease depends weakly on the type of noise 
in the probe, and for x _1 > 0.05 it does not exceed 6%. 
In this sense, the position of the peaks is not affected 
significantly by the presence of fluctuations in the probe, 
and this enables us to deduce the dipole moment for the 
transition between the Auger states |2) and |3). Indeed, 
the pump is Fourier limited and in principle one can have 

control on the peak intensity I d . Hence, in view of Eq. 

(|2"rj)) . one can deduce fi23 by estimating Q d °^ through the 
separation of the peaks in the average stochastic signal 



FIG. 9. (Color online) DR arrangement I with stochastic 
probe and Fourier-limited pump, (a-c) The separation of the 
peaks in the average stochastic signal (Q2), normalized to 
the separation in the signal for Fourier-limited probe of the 
same shape (approximately equal to 19.2F2), as a function of 
X _1 for various noises, (d-f) The depth of the doublet in the 
average stochastic signal (given by Eq. I28[) . as a function of 
X -1 for various noises. Other parameters as in Fig. [8] 



(QaC^s)}) f° r Ad = 0. As discussed in the appendix, 
for Fourier-limited probe and pump pulses, the separa- 
tion of the peaks is not precisely equal to 20,^ , but it is 
very close [~ 8% smaller for the parameters of Fig. Eta)]. 
Thus, obtaining fi23 from the splitting in (Q2(A S )) will 
unavoidably lead to an error due to the time-dependence 
of the pump and the probe, as well as an error due to the 
fluctuations in the probe. Assuming that the Fourier- 
limited profile of the pump is known, the former error 
could be eliminated only if one can measure the aver- 
age intensity profile for the probe (I s (t)} together with 
(QaC^s)}- br an y case, we believe that for judiciously 
chosen probe/pump pulse durations, and coherence time 
for the probe, this approach allows one to measure the 
dipole moment between Auger states with at most 20% 
of error. To the best of our knowledge, this is the only 
way to the measurement of the dipole moment between 
two Auger or autoionizing states, which provides a very 
sensitive test of theoretical models for states embedded 
in continua. 

Besides the splitting, the visibility of the doublet also 
depends on its depth which can be quantified by 



V 



(Q 3 -(A fl )) r 



(Qi(A.)> r 



(Oj(A.)) t 



(28) 



and is also plotted in Figs. [9f d-f ) as a function of x _1 
for various noises. For x ~ lj V" is close to 1 and de- 
creases with decreasing x _1 - The critical value of x" 1 
below which we observe a rapid decrease of the depth is 
approximately 0.1, and depends weakly on the type of 
the noise in the probe. It basically corresponds to FEL 
bandwidth 7 w Qd/2. As a general remark, note that 



12 




-(b) 




6 12 18 

X 



FIG. 10. (Color online) DR arrangement I with stochastic 
probe and Fourier-limited pump, (a) The FWHM of the 
two peaks in the average stochastic signal (Q 2 (A S )), is plot- 
ted as a function of x, f° r various types of noises. Other 
parameters as in Fig. [8] (b) As in Fig. UOf a). but for 
Gaussian correlated noise only, and various combinations of 
atomic and field parameters. (Circle) Gaussian profile for 
(I a (t)), t s — 4.5r 2 ^ 1 . r3 = r 2 ; (Square) Gaussian profile for 
(I s {t)), t s = 4.5rr\ T 3 = 0.5F 2 ; (Star) Gaussian profile for 
(I s {t)), t s = 3^, r 3 = r 2 ; (Triangle) profile 2 for (J s (t)), 
t s = 4.5r^" , r$ = r 2 . Other parameters: Gaussian pump 

pulse, r d = fii 0) 
5000 realizations. 



= o.ir 2 , n { d 0) = ior 2 , i> = 32rj\ 



there is a rather large regime of values for where the 
depth and the splitting do not vary appreciably. Hence, 
none of them seems capable of resolving reliably different 
values of x, providing thus information about the actual 
bandwidth of the probe. The question arising therefore 
is whether the width of the peaks are more sensitive to 
variations of Xi an d may thus provide more information 
in this respect. 

In Fig. [TUta) we plot the mean value of the FWHM for 
the two peaks in the average stochastic signal, as a func- 
tion of x, for various types of noises in the probe pulse. 
Clearly in all of the cases, we have a linear dependence 
on x- The corresponding slope is the same for Gaussian 
and exponentially correlated noise, whereas it is smaller 
for noise with correlations that drop as a hyperbolic se- 
cant (see tabic QI . According to Fig. HUT b) , the linear 
dependence and the slope depend only weakly on the pro- 
file of the average probe intensity (I s (t)) and the decay 
rate of state |3). However, the slope does seem to have 
an inversely proportional dependence on the duration of 
(I s (t)} (this was to be expected of course, in view of the 
discussion in the appendix). These observations suggest 
a strong correlation between the FWHM of the peaks in 
the average stochastic signal and the bandwidth of the 
probe Aw s , from the regime of Fourier-limited probes 
(x <C 1) all the way up to stationary probes 1). 



2. Arrangement II: stochastic pump 

For this arrangement, the transition |1) -f-> \2) is 
strongly driven by stochastic pulses, whereas the Fourier- 
limited probe is applied on |2) <-> |3). The splitting in- 
duced in the lower transition is monitored by keeping 
track of the Auger electrons from the decay of state |3). 
The average stochastic signal is depicted in Fig. [TTJ as a 
function of the probe detuning Ad, and for three differ- 
ent detunings of the pump. Clearly, the stochastic pump 
has a dramatic effect on the doublet, as compared to the 
case of Fourier-limited pump (dashed curves). The peaks 
of the doublet broaden significantly even for values of T c 
comparable to r s , whereas at the same time their sepa- 
ration is reduced. 

The effect of fluctuations in this arrangement are far 
more severe than in arrangement I, because they per- 
tain to the pump pulse, and not to the probe. The 
atoms experience different intensities at different times 
within a single pump pulse, and the induced splitting 
also changes with time. In the picture of eigenstates (see 
appendix), instead of two states |±) with well-defined en- 
ergies, one has essentially two manifolds of |+) and |— ) 
states, and there are no clear resonances for the probe 
to resolve. Hence, the distribution of (^(A^) exhibits a 
rather strong background and the visibility of the dou- 
blet is significantly reduced. The background becomes 
stronger as we increase x 101 fixed t s (i..e, with increas- 
ing bandwidth or decreasing T c ), leading to severe quali- 
tative and not only quantitative deviations from the cor- 
responding signal for Fourier-limited pump. By contrast 
to arrangement I, in the present arrangement the width 
of the peaks in the doublet are comparable to their sep- 
aration, which reduces even further the visibility of the 
doublet. Hence, it is rather difficult to extract from the 
doublet, accurate information e.g., about the bandwidth 
of the FEL pulses. In this respect, therefore the previous 
arrangement is far more useful for spectroscopic applica- 
tions. 



V. SUMMARY 

We have discussed the effects of field fluctuations on 
the frequency response of single and double Auger res- 
onances driven by FEL radiation. We have found that, 
in the case of single resonance the fluctuations tend to 
enhance power broadening phenomena, leading thus to 
unconventional distributions for the total yield of Auger 
electrons with respect to the detuning of the driving field 
from resonance. In the case of double resonance, when 
the FEL radiation plays the role of the probe, the total 
yield of the Auger electrons exhibits the characteristic 
Autler-Townes doublet when the bandwidth of the source 
does not exceed the Rabi frequency of the pump. The 
widths of the peaks have been shown to be strongly cor- 
related with the bandwidth of the FEL pulses, whereas 
their separation allows for the estimation of the dipolc 
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FIG. 11. (Color online) DR arrangement II with stochastic 
pump and Fourier-limited probe. The average stochastic sig- 
nal is plotted as a function of the probe detuning Ad, for 
increasing values of the ratio \, and for pump detunings (a) 
A s = 0; (b) A s = 0.ir 2 ; (c) A s = 2F 2 . The thin dashed 
curves correspond to a Fourier-limited pump (i.e., ^ < 1), 
Synchronized Gaussian profiles have been chosen for the aver- 
age pump intensity (I s (t)) as well as for the probe. Other pa- 
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FIG. 12. (Color online) DR arrangement I with Fourier- 
limited pump and probe. The Auger signal Q 2 is plotted as a 
function of the detuning A s , for two different ratios of fi d °' /r 2 
and pump detunings (a) Ad = 0; (b) Ad = 2F 2 . The dashed 
curves correspond to synchronized Gaussian probe and pump 
pulses with durations t s = i.ST^ and Td = BT^ , respec- 
tively. The solid curves correspond to constant pump and 
probe fields lasting for time equal to Tf . For the sake of com- 
parison, the signal has been normalized to its maximum value. 
Other parameters: fii 0) = 0.ir 2 , r 3 = T 2 , T f = 32^ 1 . 



moment between the Auger resonances involved in the 
double-resonance setup. Finally, when FEL radiation 
is used to pump the double-resonance, field fluctuations 
have been shown to have a dramatic effect on the ac Stark 
splitting. The modelling of fluctuations and the parame- 
ters used throughout this work follow closely real experi- 
mental situations, and thus our findings can be tested in 
practise. 

Although our discussion and quantitative illustration 
of the physics has been in the context of Auger reso- 
nances, the formalism and treatment are equally valid 
for any type of decaying resonances, in any system. In 
closing, we should mention that we are well aware of the 
effects of interaction-volume expansion in any situation 
involving strong radiation, which by necessity is focused. 
This is an instrumental effect which is apt to affect the 
observed line shapes, and as such needs to be taken into 
consideration in the interpretation of experimental data. 
It does, however, depend on the particular focusing ge- 
ometry pertaining to a given experiment, but the relevant 
theoretical tools are known [Tl(. In any case, having re- 
sults without the volume effect provides a point of cali- 
bration to be taken into consideration in the estimation 
of its expected influence. 
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Appendix: DR with Fourier-limited pulses 

The problem of DR with Fourier-limited pulses has 
been studied extensively in the past, and the reader may 
refer to related literature (e.g., see [22, l35r!4(ij and ref- 
erences therein). To facilitate the discussion of Sec. IIV1 
we summarize briefly here fundamental aspects of this 
problem. 



1. Arrangement I: probe in the lower transition 

Consider the case where the upper transition is 
strongly driven, and the lower one is weakly driven by 
the probe. The strong driving of the upper transition 
results in a splitting referred to as ac Stark splitting. In 
the case of stationary fields (i.e., for time independent 
£s(d)(t)), one can readily show this splitting by diago- 
nalizing the Hamiltonian corresponding to the strongly- 
driven transition |2) o |3). There are two eigenencrgics 
given by 

with the corresponding eigenstates |±) = =Fsin(0±) |2) + 
cos(9±) |3), where tan(0±) = z f£ld/w±. The induced ac 
Stark splitting is simply uj + — cj_ = J + 4£1^, and 
can be monitored by following the absorption of a weak 
probe, as we vary its frequency around the lower reso- 
nance |1) <H- 1 2). The DR setup in this picture involves 
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FIG. 13. (Color online) DR arrangement I with Fourier- 
limited pump and probe, (a-c) Signal Q 2 at the end of the 
sequence of probe/pump pulses shown in the upper panel, as 
a function of the detuning A s . The probe follows the pro- 
file 2, discussed in Sec. [TH and its FWHM is approximately 
7.5F2 1 . The pump pulse is taken Gaussian with duration: (a) 
Td — iY^ 1 ', (b,c) Td — GY^ 1 . For (a) and (b) the probe is de- 
layed relative to the pump by approximately G.&Y^ 1 , whereas 
for (c) the probe precedes by approximately 'YY -1 . Other pa- 
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FIG. 14. (Color online) DR arrangement II with Fourier- 
limited pump and probe. The Auger Signal Q3 is plotted 
as a function of the detuning Ad, for two different ratios of 
fi£°'/r2, and for pump detunings (a) A s = 0; (b) A s = 0.ir 2 . 
The dashed curves correspond to Gaussian probe and pump 
pulses with durations Td = 3r^~ and t s = BF^" 1 , respec- 
tively. The solid curves correspond to constant pump and 
probe fields lasting for time equal to Tf. For the sake of com- 
parison, the signal has been normalized to its maximum value. 
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two distinct absorption paths |1) — > |±) for the probe, 
with the ratio of the corresponding couplings determined 
by the overlaps i.e., fi = |(1 1-) |/|(1 |+) |. For A d = 0, 
pb = 1 and thus the two absorption channels are equiv- 
alent, whereas for > the channel |1) — > |+) dom- 
inates over the channel |1) — > |— ), and the opposite is 
true for A^ < 0. The absorption of the probe with vary- 
ing A s can be monitored directly (e.g., via a detector), 
or indirectly e.g., by examining the Auger electrons from 
level |2); larger absorption yields stronger signal. Such a 
signal in presented in Fig. I12| for two different Rabi fre- 
quencies and pump detunings (see solid lines). Clearly, 
Q2 exhibits a doublet with the two peaks located at u)± , 
and the highest peak being the one at cj + for A^ > 0. 
Moreover, in the regime of A s = the signal is negligible. 
In the spirit of EIT, this corresponds to a transparency 
window, where the absorption of the probe is reduced 
significantly to frequencies around W21 ■ 

The previous qualitative description also holds for 
pulsed pump and probe. In this case, however, the atoms 
experience different intensities at different times, and this 
results in a considerable broadening of the peaks, as well 
as a small decrease of the separation of the peaks in the 
doublet (see dashed curves in Fig. [12]). The FWHM of 
the peaks in the doublet acquires a dependence on the 
the pulse duration, in addition to the dependence on the 
lincwidths of the states involved. In contrast to station- 
ary fields, there are many additional parameters in the 
problem pertaining to the details of the pulses (such as 
shapes, durations, peak values, etc) as well as their rel- 
ative position. As a rule of thumb, in order for a clear 
doublet to be observable, the atoms have to be "dressed" 
by the pump field before the probe rises, so that the split- 
ting has been induced. The conditions under which such 
a requirement is fulfilled depends on the particular pulse 
shapes under consideration. Roughly speaking, the du- 
ration of the pump must be sufficiently longer than the 
duration of the probe pulse (i.e., > r s ), and any rela- 
tive delay between the pulses should be sufficiently small 
so that the overlap between the pulses is not reduced sig- 
nificantly. As depicted in Fig. [T3l if these conditions are 
fulfilled, a clear doublet can be observed irrespective of 
the actual profile of the probe, whereas in the opposite 
case the doublet may be distorted by additional artifacts 
(essentially the probe resolves various irrelevant details 
during the rise and the fall of the pump pulse). The 
fact that the observability of the doublet does not seem 
to depend on the form of the pulse shapes is of particu- 
lar interest for the discussion in Sec. IIV1 because when 
SASE-FEL pulses arc involved in the DR, one cannot 
predict in advance the precise profile for the average in- 
tensity (I a (t)). 



2. Arrangement II: probe in the upper transition 



In this case, the diagonalization has to be performed on 
the part of the Hamiltonian that pertains to the strongly 
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driven transition |1) — > |2). The above expressions for 
the eigenergies and eigenstates hold, with the replace- 
ments fl d -> fl s , A d -> A s , |2) -> |1) and |3) -> |2). By 
contrast to the previous arrangement, the absorption of 
the probe now involves the two distinct paths |±) — > |3), 
and can be monitored e.g., by following the Auger elec- 
trons from level |3). The couplings for the two absorption 
channels |±) — > |3) depend on the overlaps |(3 |±) |, with 
the corresponding ratio being equal to 1/fi. As a result, 
in the case of off- resonant drive (i.e., for Ad ^ 0) the 
position of the highest peaks in the asymmetric doublet 



is opposite to the one in arrangement I. However, the 
signal Q3 is substantially weaker than in the previous ar- 
rangement, because of the admixture of state |1) with the 
decaying state |2) in the eigenstates. The strong driving 
of 1 1) <-> 1 2) results in a fast depletion of the population 
by optical pumping, and thus to significant losses that do 
not contribute to Q3. For the same reason, in the case of 
off-resonant drive, the shorter peak in the doublet seems 
to be suppressed severely, even for a small increase of 
|A S |. Finally, the effects of time-dependent pump and 
probe on the ac Stark splitting and its observability are 
analogous to the ones in arrangement I. 
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